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To calculate the magnetic ground state of nanoparticles we present a self-consistent first princi- 
ples method in terms of a fully relativistic embedded cluster multiple scattering Green's function 
technique. Based on the derivatives of the band energy, a Newton-Raphson algorithm is used to 
find the ground state configuration. The method is applied to a cobalt nanocontact that turned 
out to show a cycloidal domain wall configuration between oppositely magnetized leads. We found 
that a wall of cycloidal spin-structure is about 30 meV lower in energy than the one of helical spin- 
structure. A detailed analysis revealed that the uniaxial on-site anisotropy of the central atom is 
mainly responsible to this energy difference. This high uniaxial anisotropy energy is accompanied 
by a huge enhancement and anisotropy of the orbital magnetic moment of the central atom. By 
varying the magnetic orientation at the central atom, we identified the term related to exchange 
couplings (Weiss-field term), various on-site anisotropy terms, and also those due to higher order 
spin-interactions. 



INTRODUCTION 



As magnetic storage devices approach a physical limit 
of a single atom, the investigation of nanoclusters has 
become one of the most important subjects in mag- 
netism. Recent developments in nanotechnology permit 
the construction of clusters with well-controlled struc- 
tures and enable the measurement of various magnetic 
properties at the atomic scale. Probing the Kondo res- 
onance in terms of low temperature scanning tunneling 
spectroscopy Heinrich et alX determined the spin-flip en- 
ergy of a single manganese atom on a nonmagnetic sub- 
strate, while Wahl et al£ were able to estimate the ex- 
change coupling between Co atoms on a Cu(001) surface. 
Atomic scale contacts can be fabricated by using electro- 
migrated break junctions where the size of a macroscopic 
contact between two leads can be reduced down to a sin- 
gle atom. Neel et al& studied the transition from the 
tunneling to the contact regime by moving the STM tip 
closer to the surface adatom, and an enhanced Kondo 
temperature was found. In conjunction with the Kondo 
effect, Calvo et alA found a Fano resonance for ferromag- 
netic point contacts indicating that the reduced coordi- 
nation can dramatically effect the magnetic behavior of 
nanoclusters. 

Experiments on atomic-sized contacts of ferromagnetic 
metals generated by mechanically controllable break 
junction (MCBJ) revealed magnetoresistance (MR) ef- 
fects of unprecedented size£~— There are various mech- 
anisms to this huge MR discussed in the literature: de- 
pending, e.g., on the micromagnetic order of the sam- 
ple controlled by the size of the applied field, atomically 
enhanced anisotropic MR (AAMR), giant MR (GMR), 
tunnel MR (TMR) or ballistic MR (BMR) effects can 
be established. 8 In particular, based on ab initio calcula- 
tions, the AAMR has been shown to emerge in wire-like 



transition metal nanocontacts and related to the giant 
orbital moment formed at the central atom^ 

Ab-initio calculations on magnetic nanostructures are 
useful for a clear interpretation of experimental re- 
sults and to attain better understanding of the un- 
derlying physical phenomena. Several methods to de- 
termine complex magnetic ground states of nanopar- 
ticles from first principles are based on a fully un- 
constrained local spin-density approximation (LSDA) 
implemented within the full-potential linearized aug- 
mented plane- wave (FLAPW) method 1 — or the projector 
augmented-wave (PAW) method^ Unconstrained non- 
collinear magnetic calculations are also performed within 
a tight-binding approach^ using the tight-binding lin- 
earized muffin-tin orbital (TB-LMTO) metho d 13 ' 14 or 
the Korringa-Kohn-Rostoker (KKR) method. 15 Spin- 
orbit coupling (SOC) has an important role in the for- 
mation of different magnetic states via magnctocrys- 
talline anisotropy and Dzyaloshinsky-Moriya (DM) inter- 
actionS; 1 - SOC is usually treated as perturbation or by 
directly solving the Dirac equation. The latter concept 
is applied in studies relying on ab-initio spin-dynamics 
in terms of a constrained LSDA by means of a fully rel- 
ativistic KKR methodJI"ia 

In bulk ferromagnets the formation of a domain wall 
is governed by a competition between the exchange and 
anisotropy energies^ and the typical interface between 
the magnetic domains is the Bloch wall where the mag- 
netization remains perpendicular to the axis of the wall. 
In thin films with easy plane anisotropy, a Neel wall 
is formed with atomic magnetic moments lying in the 
plane of the film, however, DM interactions can give 
rise to domain walls with out-of-plane magnetization and 
well-defined rotational sense d 1 ' 22 In a geometrically con- 
strained system the structure of a domain wall is mainly 
determined by the geometry irrespective of the exchange 
and anisotropy energies. 2 - 3 Thermal effects play an ad- 
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ditional role and can lead to new types of domain walls 
beyond the usual restriction of constant magnetization 
magnitude. 24 However, for a deeper understanding of the 
magnetic properties of nanocontacts, models based on 
first princples calulations are of pronounced importance. 

In the present work, a domain wall through a point- 
contact between (001) surfaces of fee Co is studied, where 
the magnetizations are aligned in the (110) and the (110) 
directions in the leads. It should be noted that Co ex- 
hibits a hep structure in bulk, however, as thin film it 
often displays a fcc-related geometry. We apply a fully 
rclativistic embedded cluster Green's function technique 
based on the KKR method (EC-KKR).^ 5 . Using gradi- 
ents and second derivatives of the band energy related to 
the transverse magnetization, a self-consistent Newton- 
Raphson method is developed to find the ground state 
configuration of the domain wall. An enhancement of the 
magnetic anisotropy energy has been established theoret- 
ically in atomic scale junctions even for elements that are 
nonmagnetic in bulk. 26 In agreement with this finding, 
our results reveal that the central atom with the lowest 
coordination number has the main contribution to the 
magnetic anisotropy of the contact. To highlight the re- 
lationship between the obtained cycloidal domain wall 
configuration and the magnetic anisotropy, the orienta- 
tional dependence of the band energy of the point-contact 
is analyzed in details. 



II. COMPUTATIONAL DETAILS 

Our model of the atomic-sized point contact has been 
built from Co atoms forming two identical pyramids fac- 
ing each other between (001) interfaces of fee Co as it 
is shown in Fig. [lja). The distance between the central 
atom and its neighbors was chosen identical to the fee 
nearest neighbor distance, a, of 2.506 A. Note that this 
geometrical model is the same as the one labelled by C2 
in Ref. 0, except that they studied a break-junction be- 
tween bec Fe surfaces. In order to mimic the contraction 
and expansion of the contact, the normal to plane dis- 
tances in the vicinity of the central atom have been scaled 
by a factor, hereinafter denoted by x, between 0.85 and 
1.15, see Fig. QJa). A host system assembled of two op- 
positely magnetized semi-infinite Co leads and separated 
by 7 layers of empty spheres (vacuum) is considered. The 
embedded cluster in the EC-KKR calculations consisted 
of 29 (9 + 4 + 1 + 1 + 1 + 4 + 9) Co atoms forming the con- 
tact by substituting empty spheres in the vacuum layers, 
16 + 16 Co atoms from the Co surfaces adjacent to the 
contact, and we also included 80 empty spheres in the 
vicinity of the Co atoms in the contact to let the electron 
density relax around the cluster, see Fig.QJb). 

First, the electronic structure of the host was cal- 
culated in terms of the fully relativistic screened 
KKR method applying the surface Green's function 
technique. 27,28 Then the electronic structure of the con- 
tact has been determined within the EC-KKR method, 25 
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Figure 1. (Color online) (a) The geometry of the contact 
viewed from the (110) direction. The leads are depicted as 
dark (blue) rectangles, the cobalt atoms forming the contact 
are represented by gray (orange) circles, a denotes the near- 
est neighbor distance in the fee structure. The length of the 
contact is tuned via x = 0.85, 0.90, 0.95, 1.00, 1.05, 1.10, 
and 1.15. Note that only the marked distances were scaled, 
(b) Sketch of the embedded cluster. Dark (blue) circles: se- 
lected atoms of the cobalt leads, gray (orange) circles: cobalt 
atoms in the nanocontact, and empty circles: empty spheres 
around the contact. The directions of magnetization in the 
leads are marked by dark (blue) arrows. 



in which the scattering path operator (SPO), correspond- 
ing to a finite cluster, C, embedded into a host system can 
be obtained from the following equation, 



r c ( £ ) = (tc 1 ( £ )-t h 1 (e)+r h - 1 ( e ))' 



(1) 



where th(e) and T\(e) denote the single-site scattering 
matrix and the SPO matrix for the host confined to the 
sites in C, respectively, while tc denotes the single-site 
scattering matrices of the embedded atoms. The calcu- 
lations for both the host and the cluster were performed 
within the local spin-density approximation (LSDA), 29 
by using the atomic sphere approximation (ASA) and 
f ma i = 2 for the angular momentum expansion. 

A fully unconstrained extension of the relativistic EC- 
KKR method is used to find the magnetic configuration 
of the point contact. The evolution of the atomic mag- 
netic moments is treated in a semi-classical manner simi- 
lar to molecular dynamics, whereby, in spirit of the mag- 
netic force theorem,— the driving force is calculated as 
the derivative of the band energy, 



(e — £f) n{e) de 



N(e)6e, (2) 



with respect to the transverse change of the exchange 
field, where £f is the Fermi energy, while n(e) and N(s) 
stand for the density of states (DOS) and for the in- 
tegrated DOS, respectively. In the multiple scattering 
formalism the exchange field enters the electronic struc- 
ture via the single-site scattering matrix, ti. The first 
and higher order changes of the ti matrices as well as the 
derivatives of the band energy can straightforwardly be 
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calculated in the local frame of reference introduced at 
all sites of the cluster, where the direction vector cr^ of 
the magnetization at site i, and the two transverse vec- 
tors, e^i and e^, form a right handed coordinate system 
as shown in Fig. [2j The first and second order change of 
the single site scattering matrix at site i with respect to 
rotations by Acf>i a around the transverse axes ei a can be 
given by the following commutator formulas, 



At) 1 ' = i[e ia J,ti}A(t)i a , 

(2) 

At] = -[e ia J, [e^J,fi]]A0 iQ A^, 



(3) 
(4) 



where J is the matrix representation of the total angular 
momentum operator and a, (3 £ {1,2}. Following Ref . l3lL 
the first and second derivatives of the band energy can 
then be expressed as 

d 6F 

—A = - Re / Tr {t u [e ia J, m»]} de, (5) 

— oo 

£F 

d^E If 
Q , J", = Im / Tr {r^ [e j/3 J, mj]Tji[e ia J, rrn}} de 

0<PiaO(Pjp 7T J 

— oo 

+ ^-Im / Tr {Tu[e ia J, [e,^J,mi]]} de, (6) 



where m ; = t~ and r.y is the block of the SPO matrix 
between sites i and j. Note that for brevity we dropped 
the energy arguments of the corresponding matrices in 
Eqs. (j3H6j) . In the spirit of a gradient minimization, ro- 
tating the exchange field by a small amount around the 
torque vector at each sites, 



dE h 



e i2 - 



dE h 



l>,2 



(7) 



the magnetic configuration gets closer to the local mini- 
mum of the energy, however, the convergence is very slow. 
In order to speed up this procedure, a Newton-Raphson 
iteration scheme has been applied, where the inverse of 
the second derivative tensor, also referred to as the Hes- 
sian, Eq. is used to estimate the angle of rotations 
around the torque vector given by Eq. ([7]). The eigen- 
values of the Hessian also provide information about the 
stability of the a configuration with zero torque: if the 
Hessian is a positive or negative definite matrix then the 
given configuration is stable or unstable state of equi- 
librium, respectively. Once the Newton-Raphson itera- 
tion has converged, new effective potentials and exchange 
fields are generated and the procedure is repeated un- 
til the effective potential converged and the torque in 
Eq. (J7J is decreased below a predefined value of, typi- 
cally, 10~ 4 mcV. 

The starting magnetic configuration for the above op- 
timization procedure has been determined by Monte 
Carlo simulated annealing based on a simple isotropic 




Heisenberg model, Ti. 



|X)iyy Jij (T i cr j> where J lJ is 



Figure 2. (Color online) Sketch of the local frame of reference: 
the unit vector (Ti is parallel to the magnetization at site i, 
while the unit vectors en and ei2 point into the transverse 
directions. Rotations around these axes by (j>n and 4>%2 are 
also indicated. 



the isotropic exchange coupling between sites i and j. 
The coupling coefficients between the atomic moments 
were calculated by using the torque method proposed by 
Liechtenstein et al£^ The exchange couplings have been 
calculated in a ferromagnetic spin-configuration parallel 
to the (100) direction. 

In order to avoid the difficulties arising from the con- 
tinuous degeneracy of the spin-states in a Heisenberg 
model, the magnetization on the central atom was fixed 
normal to the bulk magnetization. Considering the in- 
version symmetry of the point-contact, only the (110) 
and the (001) directions are consistent with the (con- 
strained) magnetic ground-state of the system. In the 
first case, the magnetic moments at all sites (layers) re- 
main within the (001) plane, i.e., normal to the axis of 
the point-contact, therefore, in the following this spin- 
configuration will be termed as a helical domain wall. In 
the second case, all the spin moments are confined to 
the (110) plane, thus, we shall call this case the cycloidal 
domain wall. Note that the helical and cycloidal spin- 
configurations closely resemble the Bloch and Neel types 
of domain walls well-known in bulk and thin-film mag- 
nets, respectively. Since, however, these types of domain 
walls are distinct through the magnetostatic energy, to 
avoid confusion we skipped using the traditional termi- 
nology. 



III. RESULTS AND DISCUSSIONS 

A. Domain wall configurations 

Self-consistent potentials and exchange fields have 
been first determined for both the cycloidal and the he- 
lical domain walls and the Newton-Raphson iterations 
were started from both initial configurations. Interest- 
ingly, when starting from a helical spin-configuration, 
the gradients, Eq. (O, were initially zero, but the Hes- 
sian had a negative eigenvalue indicating that the heli- 
cal spin-configuration belonged to a saddle point of the 
energy surface. Throwing the system off this saddle 
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Figure 3. (Color online) The cycloidal spin-configuration ob- 
tained for the unstretched contact (x = 1). The lengths of 
the arrows, indicated also with color coding, are proportional 
to the size of the spin magnetic moments. 



point, the Newton- Raphson iterations converged to the 
cycloidal spin-configuration. Thus, independent of the 
starting configuration, the magnetic state of the nano- 
j unction converged to the cycloidal wall structure for the 
stretching range considered. In Fig. |3] the ground-state 
cycloidal wall configuration is displayed for x = 1. 

At sites within the same geometrical layer, we ob- 
tained fairly similar orientations for the magnetic mo- 
ments, therefore, the shape of the domain wall can well 
be characterized by orientations determined as an av- 
erage within layers. In Fig. |4] such a profile is shown 
for x = 1 in terms of polar angles, Remarkably, 
the well-known analytical form, i3(z) = — § tanh(2z/d w ) 
could be well fitted defining, thus, the width of the do- 
main wall, d w . This fit is also shown in Fig. 21 We note 
that following Ref. [13 the analytical form of a constrained 
wall profile should be better described by Jacobian sine 
functions. However, testing this alternative approach re- 
sulted in to a relative deviation of less than 0.5 % in the 
fitted domain wall thicknesses. 

The change of the width of the domain walls against 
the length of the point contact is shown in Fig. [5] For a 
clear interpretation, the width of the walls is normalized 
to the width of the domain wall for x = 1. As obvious 
from this figure, d w (x) ~ a;d w (1.00) demonstrating that 
the width of the domain walls follows the length of the 
point contact. In case of Fe2oNigo thin films it has been 
experimentally found that the constrained geometry can 
reduce the width of the Neel wall<^ The effect is more 
pronounced in ultrathin films of few atomic layers where 
the width of the domain wall can be as small as few 
nanometers in the vicinity of a step edge^ The effect 
of the reduced dimensionality is even more obvious in 
the case of a point contact. Since the exchange energy 
gain for the few atoms of the contact is small compared 
to the increase of the exchange energy of the leads, the 
domain wall can not penetrate into the substrates and the 
wall is confined to the contact. The same conclusion has 




z, distance from the central atom 

Figure 4. Polar angles averaged within a layer of cobalt 
atoms in the contact with x = 1 as the function of the 
distance from the central atom (in units of the fee near- 
est neighbor distance, a). The solid curve displays the fit, 
■d(z) = -J tanh(2z/d„). 




0.85 0.90 0.95 1.00 1.05 1.10 1.15 
x, stretching factor 



Figure 5. Width of the domain walls through the point con- 
tact as a function of the stretching factor, x. Note that 
dw(l-00) = 2.34 a, where a is the fee nearest neighbor dis- 
tance. The solid line stands for the identity function. 

been drawn by Bruno^ based on a theoretical study of a 
continuous model of domain walls in a confined geometry. 



B. Magnetic moments 

The low coordination in thin films and in nanostruc- 
tures is often accompanied by the enhancement of the 
atomic spin and orbital moments. In Fig. [5] the calcu- 
lated values of the local spin and orbital moments are 
given in a point contact with cycloidal wall configuration 
and stretching factor, x = 1. Since the orbital moment is 
found almost parallel to the spin-moment at each site, we 
presented the projection of the orbital moment to the lo- 
cal spin quantization axis. Since the contact has a mirror 
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Figure 6. Calculated atomic magnetic moments ((J,b) in 
half of the nanocontact for the stretching factor, x = 1. In 
the upper and lower panels shown are the spin and orbital 
moments, /i sp in and fiorb, respectively. For comparison, the 
spin-moments at the Co surface and in the bulk are 1.82 /xb 
and 1.67 /iB/atom, while the corresponding values of the or- 
bital moments are 0.14 /xb and 0.08 /j,b- 



symmetry with respect to the horizontal plane including 
the central atom, therefore, the moments in only one half 
of the contact are displayed. Our data fit nicely to the 
observation reported in Refs. l35l andl36l that the spin and 
orbital moments at sites with lower coordination number 
are larger then at sites with larger coordination number. 
This is, in particular, true for the central atom with coor- 
dination number of only two where the values of the spin 
and orbital moments are even larger than those obtained 
for small clusters on Pt(lll) and Au(lll) surfaces^ - — 

Fig. [7] shows the spin and orbital moments of the cen- 
tral atom as a function of the stretching ratio x, for both 
the cycloidal and the helical spin-configurations in the 
point-contact. Clearly, the spin moments are fairly insen- 
sitive to the domain wall configuration: this can easily 
be understood as the relative spin-directions are nearly 
the same in the two types of domain walls. Also, there 
is only a moderate change of the spin moment in the 
range of 2.35 < Mspin < 2.49 /j,b for the stretching ra- 
tios under consideration. These values compare well to 
Mspin = 2.15 /ib arid fJ, S p ln — 2.26 /is calculated for a sin- 
gle Co adatom on Pt and Au(lll) surfaces in Refs. HH 
and[H, respectively. 

The dependence of the orbital moment of the central 
atom on the stretching is more pronounced than that 
of the spin moment: in case of a cycloidal and a heli- 
cal wall it increases from about 1 /is to 2 /is and from 
0.3 /is to 1.5 fig, respectively. Similar high values of /^ rb 
for the central atom of a wire-like Fe point-contact were 
reported in Ref. H and attributed to localized atomic- 
like electronic states treated within a full Hartrcc-Fock 
scheme. It should be mentioned that for a more reliable 
description of highly localized states, the plain LSDA we 
used in our calculations should be extended with, e.g., 
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Figure 7. The spin- and orbital moments of the central atom 
as a function of the stretching. Spin moments are displayed by 
open symbols, orbital moments are displayed by filled symbols 
as calculated in the cycloidal wall (CW, squares) and in the 
helical wall (HW, triangles) configurations. 



the local self-interaction correction, LSDA+SIC^S or the 
dynamical mean field theory, LSDA+DMFT^ 

Apparently, the orbital moment of the central atom 
is systematically larger in a cycloidal wall than in a he- 
lical wall. This can be understood since these orbital 
moments correspond to different directions: in case of a 
cycloidal wall it points along the (001) directions, while, 
for a helical wall, along the (110) direction. Such a huge 
anisotropy of the orbital moment at the central atom has 
also been observed in Ref. d. According to Bruno's the- 
ory^ this large orbital momentum anisotropy is related 
to a large magnetic anisotropy energy featuring the (001) 
direction as easy axis, which clearly corroborates our re- 
sult for the preference of a cycloidal domain wall. 



C. Rotational energy of the domain wall 

The cycloidal and helical spin-configurations of the 
point contact can be transformed into each other in term 
of a simultaneous rotation of the spin-directions around 
the axis parallel to the magnetization of the leads. The 
energy along the path of this global rotation, termed as 
the rotational energy of the domain wall, was calculated 
using the magnetic force theorem, namely, from the band 
energy of the system by rotating the orientation of the 
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Figure 8. The band energy of the nanocontact with x = 1.00 
while rotating the exchange field at each atomic sites simul- 
taneously around the (110) axis. By rotating all the spins by 
90° the system goes over from the cycloidal wall (CW) into 
the helical wall (HW). The dashed line denotes the leading 
Fourier component of the band energy, —15.2 [meV] cos(2#), 
see Eq. (|8]). Note that we shifted the zero level of the energy 
to the constant term, Ko- 

exchange field at each atomic site around the (110) axis 
and keeping frozen the effective potentials and fields as 
obtained for the ground state cycloidal wall configura- 
tion. For the case of the unstretched configuration the 
results are plotted in Fig. [5] The two minima and max- 
ima of the band energy belong to the two-fold degener- 
ate cycloidal and helical domain wall configurations. The 
height of the energy barrier between the two ground state 
cycloidal spin-configurations is 32.0 meV. Similar behav- 
ior has been found for the whole stretching range of the 
point contact. The energy differences between the two 
types of domain wall as a function of the stretching ratio 
are displayed by diamonds in Fig. [3J 

Due to time reversal symmetry, the magnetic 
anisotropy energy has a periodicity of tt, but it does not 
comply with a usual cos 2 (#) dependence. To explore this 
deviation we performed the Fourier expansion, 

oo 

E h (d)=K + ]T K k cos(k9), (8) 

fc=2,4,... 

for the contacts with different stretching. Note that be- 
cause of the inversion symmetry of the contact E-^(9) = 
Eh(TT — 0) applies, therefore, the sm(k9) (k — 2,4, . . .) 
terms do not appear in the expansion, Eq. ([5]). We sum- 
marized the Fourier coefficients, Kk, in Table|H We found 
that in each case the term K% cos(2#) adds the largest 
weight to the rotational energy of the domain wall. The 
next term K4COs(49) is quite significant for x > 0.95, 
but it drops for smaller stretching. Interestingly, in the 
stretching range of x < 0.90 the fc = 6 term overweights 
the one of k = 4, whereas in the complementary range 
the k = 6 term is negligible. It should be mentioned that 



Table I. The k = 2, 4 and 6 Fourier coefficients (in units of 
meV) of the rotational energy of the point-contact, Eq. (JS|, 
as a function of the stretching parameter, x. 



x_ k = 2 k = 4 k = 6 

0.85 -6.3 0.15 0.397 

0.90 -10.0 0.36 0.499 

0.95 -13.6 1.40 0.298 

1.00 -15.2 2.17 -0.040 

1.05 -15.1 2.35 -0.122 

1.10 -14.4 2.24 -0.083 

1.15 -13.2 1.92 0.025 



the k > 8 terms of the Fourier expansion have practically 
vanishing weight. 

D. Magnetic anisotropy of the central atom 

As we have seen in Sec. MI Bl the central atom of the 
contact exhibits a huge orbital moment anisotropy that 
should be accompanied by a large magnetic anisotropy 
energy. For that reason, we analyze the band en- 
ergy of the point-contact, E\,((t), with a denoting the 
spin-orientation at the central atom, whereas the spin- 
orientations of all the other sites in the contact are kept 
fixed as obtained in the ground-state cycloidal wall con- 
figuration. 

Our analysis is based on an expansion of Eb(<r) in 
terms of (real) spherical harmonics, i?™(<x), 

E h (a) = J2K7 l RT('r), (9) 

with the angular momentum indices, £ = 0, 1, 2, . . . and 
— £ < m < I. Similar to the rotational energy of the do- 
main wall, we used the magnetic force theorem to evalu- 
ate Eb(cr), but here we employed Lloyd's formula^ since 
it accurately accounts for the change of the band energy 
of the whole point-contact with respect to the change 
of the spin-orientation at the central site. For the ex- 
pansion, the integration over <r was performed using a 
51 points Gaussian quadrature along the z-direction and 
a uniform mesh of 100 points in the azimuth angle, re- 
sulting in a spherical grid of 5100 points. The obtained 
coefficients are summarized in Table [TT] up to £ = 4 and 
for all the stretching ratios under consideration. Only 
the non-vanishing coefficients are presented, for clarity, 
together with the definition of the corresponding spheri- 
cal harmonics, i?™(<x). 

The absence of certain spherical harmonics in expan- 
sion Eq. can be discussed based on group-theoretical 
arguments. The function E^(cr) should be invariant un- 
der symmetry transformations, g, of the point-contact, 
Eb(c) = Eb(gcr) : including the symmetry of both the 
lattice and the given (cycloidal) spin-configuration. Re- 
garding that the spin-vectors transform as axial vectors, 
the only allowed transformation is the reflection onto the 
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Table II. Expansion coefficients K™ (in units of meV) of the band energy of the contact, see Eq. ©, according to real spherical 
harmonics R™ up to £ — 4. 



1 


m 




pwi 


x = 0.85 


x = 0.90 


x = 0.95 


x = 1.00 


x = 1.05 


x = 1.10 


x = 1.15 


1 









-240 


-247 


-235 


-212 


-192 


-176 


-159 


2 









-25.3 


-30.0 


-33.2 


-32.4 


-30.9 


-28.4 


-25.6 


2 


2 




Ia/¥(* 2 ^ 2 ) 


4.30 


2.54 


1.39 


0.51 


-0.29 


-0.92 


-1.36 


3 







(5z 3 -3 2 ) 


4.12 


3.06 


1.63 


0.71 


-0.28 


-1.43 


-2.67 


3 


2 






-0.199 


-0.093 


0.004 


0.108 


0.196 


0.267 


0.293 


4 





3 
16 


xjl (35z 4 - 30z 2 + 3) 


-0.63 


1.72 


4.60 


4.94 


5.05 


4.85 


4.32 


4 


2 


3 

8 \/ 


/f(* 2 -y 2 ) (7^ 2 -i) 


0.033 


0.125 


0.184 


0.108 


0.051 


0.001 


-0.052 


4 


4 


3 

16 1 


/f ( K 4-6xV + S/ 4 ) 


-0.007 


-0.005 


-0.018 


-0.041 


-0.088 


-0.187 


-0.345 



(001) plane: (x,y,z) — > (—x,—y,z). Thus we conclude 
that only those function can enter the expansion of E\, (<r) 
that contain even powers of the variables x and y. As seen 
from Table HIl this is fully confirmed by our calculations. 
Apparently, the expansion Eq. ([SJ shows a satisfactory 
convergence as the coefficients rapidly decrease with in- 
creasing £. An obvious exception can, however, be seen 
for K\ that for x > 0.95 overweig hts K$. Noticeably, 
among the terms with a given £, the one associated with 
the z component of the magnetization (m = 0), i.e., ex- 
cluding in-plane anisotropy, has the largest weight. 

In order to connect the above results to the rotational 
energy of the domain wall discussed in Sec. IIII CI we 
relate expansion Eq. (|9]) to a classical spin-model. Ac- 
cording to a Heisenberg model extended by relativistic 
correction o 31 ' 42 the energy in Eq. © can be expressed as 

E{tr) = £ anis ((r) + <r J cj°j - (10) 

3 

where J CJ - denote the exchange coupling tensor between 
the central site and the other sites of the contact with 
classical spin- vectors <Tj and -E an i s (<r) stands for the on- 
site anisotropy energy that, due to the tetragonal (D^h) 
point-group symmetry of the point-contact, can be ex- 
panded up to I = 4 as 

E anis (a) = Jf^H + KlRlia) + K\r\{ct) . (11) 

It is clear that the (£,m) — (1,0) term in Eq. ([9]) is 
uniquely related to the exchange coupling and, due to the 
presence of a cycloidal wall, it represents a strong Weiss 
field that orients the magnetic moment at the central 
site along the z direction. Because of the increasing dis- 
tances between the central site and the other sites of the 
contact, it is also easy to understand why this term sig- 
nificantly decreases with increasing stretching ratio. On 
the other hand, there is no (£, m) = (1, 0) term in the ro- 
tational energy of the domain wall, Eq. (JSJ), since in that 
case the relative orientation of the spins are unchanged. 
With other words, repeating the expansion Eq. © in the 



presence of a helical wall, the leading term correspond to 
the spherical harmonics oc x, with practically the same 
coefficients as listed in Table [TTI for [I, m) = (1,0). 

In relation to Eq. (jlip. the terms proportional to R^, 
i?4 and R\ in Eq. © can mainly be attributed to on-site 
anisotropy contributions to the spin-Hamiltonian, how- 
ever, the effect of higher order spin-interactions can not 
be ruled out. The second-order uniaxial anisotropy coef- 
ficients, if ° > are negative in the whole range of stretching, 
favoring thus a normal-to-plane direction. Remarkably, 
the magnitude of is around 30meV, with a maxi- 
mum of \K®\ = 33.2 meV at x = 0.95. This value should 
be compared to some results communicated in the lit- 
erature: Etz et al^ and Bornemann et calculated 
5.3 meV and 4.76 meV, respectively, for the MAE of a Co 
ad-atom on Pt(lll) surface, while, including orbital po- 
larization, Gambardella et al. 4 - obtained 18.45 me V for 
the same system. In a similar geometrical confinement 
of an atomic scale junction, W and Ir turned out to be 
magnetic with a magnetic anisotropy energy comparable 
to our values ^ 

From Fig.[S]and Table|T]we inferred that the rotational 
energy of the domain wall is dominated by the uniaxial 
magnetic anisotropy term proportional to cos 2 9 = z 2 . 
In Fig. [5] the energy differences obtained between the he- 
lical wall configuration and the ground state cycloidal 
wall configuration is plotted as a function of the stretch- 
ing factor, together with that provided by the uniaxial 

anisotropy of the central atom, |y^i^2- The values of 

AE from the two calculations agree well for x > 0.95, 
while for more squeezed contacts the uniaxial anisotropy 
of the central atom overestimates the energy difference 
between the different types of domain walls. Neverthe- 
less, we can in general conclude that the main driving 
force of the formation of a cycloidal domain wall is a gi- 
ant uniaxial on-site magnetic anisotropy at the central 
atom: in the cycloidal wall the magnetic moment of the 
central atom is parallel to the easy axis, while in the 
helical wall configuration it lies within the hard plane. 
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Figure 9. Diamonds: Calculated energy differences between 
the helical and cycloidal domain walls, Bhw — Ecw, circles: 
on-site uniaxial magnetic anisotropy energy of the central 
atom (see text) as a function of the stretching parameter, x. 
Thin lines serve as a guide for the eye. 



Finally, we briefly comment on the terms correspond- 
ing to (£,m) = (2,2), (3,0), (3,2) and (4,2) in Ta- 
ble [TTJ Since these terms are not invariant under trans- 
formations of the D^h point-group, they can not be ac- 
counted for the on-site anisotropy terms. In terms of a 
spin-model, these terms should, therefore, be related to 
higher order spin-interactions. The (t,m) = (2,2) term 
can, e.g., be identified as the consequence of biquadratic 
interactions, 4 - J2i B C i(cr<Ti) 2 , while the £ = 3 terms to 
triquadratic interactions, 47 JN T C i(er(Ti) 3 . Four-spin in- 
teractions have been explicitely calculated and proved to 
give significant contributions to a spin-Hamiltonian of Cr 
trimers deposited on Au(lll) surface by Antal et al.^ 
but recently their presence was highlighted even in bulk 
magnets 42 



IV. SUMMARY 

In case of deposited magnetic nanostructures the 
point-group symmetry of the system might consider- 
ably be reduced, therefore, complex magnetic states oc- 
cur naturally. Detecting and investigating such mag- 
netic states pose a challenge for ab initio calculations. 
We have developed a computational technique based 
on a self-consistent embedded cluster Korringa-Kohn- 
Rostokcr method suitable to find non-collinear ground- 
states of finite magnetic clusters. The method is ap- 
plied to determine the structure of a domain wall formed 
through an atomic scale nanocontact between two an- 
tiparallelly magnetized cobalt leads. The obtained 
ground state is a cycloidal domain wall which remains 
stable against squeezing or stretching the contact along 
the normal-to-plane direction. A huge enhancement, as 
well as, anisotropy of the orbital moment are found at 
the central site of the contact. The energy of the domain 
walls was explored in terms of the magnetic force theo- 
rem. Our main observation is that the formation of the 
cycloidal wall against a helical wall is primarily driven 
by the uniaxial on-site anisotropy at the central site. 
We also found effects of higher order spin-interactions as 
terms in the expansion of the band energy not complying 
with the point-group symmetry of the point contact. 
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